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ABSTRACT 
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Context. Being dominated by non-thermal emission from aligned relativistic jets, blazars allow us to elucidate the physics of extra- 
galactic jets, and, ultimately, how the energy is extracted from the central black hole in radio-loud active galactic nuclei. 
Aims. Crucial information is provided by broad-band spectral energy distributions (SEDs), their trends with luminosity and correlated 
multi-frequency variability. With this study we plan to obtain a database of contemporaneous radio-to-optical spectra of a sample of 
[ blazars, which are and will be observed by current and future high-energy satellites. 

Methods. Since December 2004 we are performing a monthly multi-frequency radio monitoring of a sample of 35 blazars at the 
antennas in Medicina and Noto. Contemporaneous near-IR and optical observations for all our observing epochs are organised, 
f Results. Until June 2006 about 4000 radio measurements and 5500 near-IR and optical measurements were obtained. Most of the 

sources show significant variability in all observing bands. Here we present the multi-frequency data acquired during the first eighteen 
months of the project, and construct the SEDs for the best-sampled sources. 
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1 . Introduction presents a large number of poorly understood details which, in 

turn, lead to a wide variety of models and call for long-term and 

H ; Blazars form a sub-group of radio-loud active galactic nuclei mu i ti . wave length campaigns capable of providing the necessary 

^ (AGN) showing extreme variability at all wavelengths, high de- observational constraints 
grees of linear polarization, strong y-ray emission, apparent su- 

perluminal motion of radio jet components, and bri ghtness tern - From the time scales of variability key information on the 

peratures exceeding the Compton limit (see e.g. [Un3[l999). AGN structure can be derived, down to linear scales or flux 

They include BLLac objects as well as quasars with flat ra- density levels not accessible even with interferometric imag- 

dio spectrum and violent variability in the optical band. The hig, but the mechanisms which are responsible for variability 

large amount of work done in the last decades led to the com- are not well understood yet. Different mod els are discussed 

monly accepted scenario in which a rotating supermassive black in the literature such as shock-in-jets (e.g. |Marscher & Gearj 

hole surrounded by an accretion disk with an intense plasma jet CUSS lAller et al.lll985HMarschenll ^96b or colliding relativis - 

closely aligned to the line of sight is responsible for the blazar tic plasma shells (e.g. |Spada et al.| |200lt |Guetta et al.| |2004|)- 

emission. Relativistic electrons produce the soft photons through Further, in the case of processing binary bla ck-hole systems (e.g. 

synchrotron emission, while hard photons are likely produced by iBegelman etal. 1980; SillanpaaeLalJ[l98 8j), rota ting helical jets 

inverse-Compton scattering. This overall scenario, however, still (e.g. |Villata et al.| |1998|; I Villata & Raiteril 1 1 999t) or helical tra- 

jector ies of plasma elements (e.g. lCamenzind & Krockenbergeil 

Send offprint requests to: U. Bach, e-mail: |1992|), the models suggest that changes in the direction of for- 

ubach@mpifr-bonn.mpg.de ward beaming introduces flares due to the lighthouse effect. 

* Current address: Max-Planck-Institut fur Radioastronomie, Auf Thus, variability furnishes important clues into size, structure, 

dem Hiigel 69, 53121 Bonn, Germany physics and dynamics of the radiating source region. A tool to 



study the radiating particles is provided by the analysis of the 
synchrotron continuum over a range of frequencies as broad as 
possible. The shape of the synchrotron spectrum gives direct 
insight into the shape of the electron energy distribution, thus 
also constraining the emission by the inverse-Compton process 
at other wavelengths. 

The most extensive blazar monitoring campaigns have been 
carried out at radio frequencies. The University of Michigan 
monitoring program at 4.8, 8.0 and 14.5 GHz has obtained data 
on over 200 sources for over three decades dAller et al.ll2003l) 
and is still continued. At higher frequencies, the Metsahovi 
Radio Observatory has reported observations at 22, 37 and 
87 GHz of 157 extragalactic radio sources ( Terasrant a et al.l 
1 998. 120041 120051) . many of which have been monitored for over 
20 years, but the m onitoring is now c ontinued only at 37 GHz. 
The Bologna group dBondiet all 19961) has observed at 408 MHz 
125 radio sources from 1975 to 1990. A sample of 550 sources 
is monitored at six radio frequencies from 1 .4 to 31 GHz quasi 
simultaneously since many years at the RATAN-600 telescope 
dKovalev et alj|2002l) . 

However, aside from some organised broad multi- 
wavelength campaigns, most of our previous studies of the 
spectral energy distributions (SEDs) of blazars were performed 
by combining measurements carried out non-simultaneously 
(over a period of one or more months), using the above men- 
tioned radio surveys and dat a from various opt i cal, near-IR , 
and X-ray observat ories (e.g. iRaiteri et all 1 19991 l200ll 120031: 
IViUata et al.l 12004 j) . Since blazars are highly variable, it is 
therewith necessary to extrapolate between non simultaneous 
points, which is not an ideal procedure. Hence, it is important 
to take simultaneous broad-band spectra that are not affected by 
the variability. Therefore, we started a monthly multi-frequency 
radio monitoring of a sample of 35 y-ray loud blazars at the 
antennas in Medicina and Noto, and organised contemporane- 
ous near-IR and optical observations for most of the observing 
epochs. This provides a valuable database of time-dependent 
SEDs, which can be used to test different jet emission models. 
The monitoring is supposed to support also the next Italian 
high-energy space missions AGILE, which will observe many 
of our monitored sources. 

In this article we will present the light curves and SEDs 
of those sources which have been monitored for at least one 
year by the end of June 2006. Further publications with a 
more detailed timing analysis and modelling of the SEDs are in 
preparation. Throughout this paper we will assume a flat uni- 
verse model with the following parameters: Hubble constant 
Hq = 71kms~' Mpc~', a pressure-less matt er content Q ? = 0-3 , 
and a cosmological constant = 0.7 (Spergel elld] 120031) . 
Cosmology-dependent values quoted from other authors are 
scaled to these parameters. 

The paper is organised as follows: Section [2] introduces the 
source sample. In Sect. [3] we will describe the data acquisition 
and reduction procedures and in Sect. |4]we will present the re- 
sulting light curves and SEDs. Peculiar variability behaviours or 
special treatments during the data analysis for single sources are 
described in Sect|5] Conclusions are drawn in Sect. [6] 



2. The source sample 

At the beginning of our monitoring we started with a sample 
of y-ray loud sources which had already been monitored in the 
optical band at the Perugia and Torino observatories since 1994 
and which were relatively bright also in the radio band (Si .4ghz > 



Uy). In view of the forthcoming launch of the high-energy satel- 
lites AGILE and GLAST, we subsequentl y added some more 
sources from the third EGRET catalog (Hart man et al.l fl999) 
with declination > -15°, which are also bright at optical (B < 
19) and radio wavelengths (Si.4ghz > 1 Jy)- The 33 sources of 
the present sample are listed in Table [TJ In this paper we present 
the data of those sources which were monitored for at least one 
year (marked in bold face in the table). The remaining ones will 
be reported in a later paper (Bach et al. in prep.). 



Table 1. The complete source list. Given are IAU name, position, 
type (Quasar, BLLac, or Seyfert) and redshift . The sources pre- 
sented in this paper are marked in bold face, the others will be 
presented in a later paper. 



IAU 


Other 


R.A. & Dec. (J2000) 


Type 


z 


Name 


Name 


[h:m:s] [d:m:s] 






0219+428 


3C 66A 


02:22:39.6 +43:02:08 


Q 


0.44 


0235+164 


AO 0235+ 16 


02:38:38.9 +16:36:59 


Q 


0.94 


0336-019 


CTA026 


03:39:30.9 -01:46:36 


Q 


0.85 


0420-014 


PKS 0420-01 


04:23:15.8 -01:20:33 


Q 


0.91 


0440-003 


NRAO 190 


04:42:38.6 -00:17:43 


Q 


0.84 


0528+134 


PKS 0528+134 


05:30:56.4 +13:31:55 


Q 


2.06 


0716+714 


S5 0716+71 


07:21:53.4 +71:20:36 


BL 


>0.3 


0735+178 


PKS 0735+17 


07:38:07.3 +17:42:19 


BL 


0.42 


0736+017 


PKS 0736+01 


07:39:18.0 +01:37:05 


Q 


0.19 


0827+243 


OJ248 


08:30:52.1 +24:11:00 


Q 


0.94 


0829+046 


OJ49 


08:31:48.9 +04:29:38 


BL 


0.17 


0836+710 


S5 0836+71 


08:41:24.3 +70:53:42 


Q 


2.17 


0851+202 


OJ287 


08:54:48.8 +20:06:31 


BL 


0.31 


0954+658 


S4 0954+65 


09:58:47.2 +65:33:55 


BL 


0.37 


1156+295 


4C +29.45 


11:59:31.8 +29:14:44 


Q 


0.73 


1219+285 


WCom 


12:21:31.7 +28:13:58 


BL 


0.10 


1226+023 


3C273 


12:29:06.7 +02:03:09 


Q 


0.16 


1253-055 


3C 279 


12:56:11.1 -05:47:21 


BL 


0.54 


1334-127 


PKS 1335-127 


13:37:39.8 -12:57:25 


BL 


0.54 


1354+195 


4C + 19.44 


13:57:04.4+19:19:07 


Sy 


0.72 


1510-089 


PKS 1510-08 


15:12:50.5 -09:06:00 


Q 


0.36 


1606+105 


4C + 10.45 


16:08:46.1 +10:29:08 


Q 


1.23 


1611+343 


DA 406 


16:13:41.0 +34:12:48 


Q 


1.40 


1633+382 


S4 1633+38 


16:35:15.4 +38:08:04 


Q 


1.81 


1641+399 


3C345 


16:42:58.8 +39:48:37 


Q 


0.59 


1730-130 


NRAO 530 


17:33:02.6 -13:04:49 


Q 


0.90 


1741-038 


PKS 1741-03 


17:43:58.8 -03:50:05 


Q 


1.05 


1807+698 


3C371 


18:06:50.7 +69:49:28 


BL 


0.05 


2200+420 


BLLac 


22:02:43.2 +42:16:40 


BL 


0.07 


2223-052 


3C446 


22:25:47.2 -04:57:01 


BL 


1.40 


2230+114 


CTA 102 


22:32:36.4 +11:43:50 


Q 


1.04 


2251+158 


3C 454.3 


22:53:57.7 +16:08:54 


Q 


0.86 


2344+092 


4C 09.74 


23:46:36.8 +09:30:45 


Q 


0.68 



3. Observations and data reduction 

3.1. Radio data 

Radio measurements were taken at the 32 m radio antennas of 
the INAF-Istituto di Radioastronomia in Medicina and Noto on 
a monthly basis at frequencies of 5, 8, and 22 GHz. At the be- 
ginning all frequencies (5, 8, and 22 GHz) were taken at both 
telescopes, but because of the better performance of the Noto 
antenna at high frequencies, the observations were split to 5 
and 8 GHz at Medicina and 22 GHz at Noto after the first 6 
months. Unfortunately, the antenna in Noto was out of opera- 
tion for the second half of 2005, because of a technical failure. 



During this time observations at 22 GHz were occasionally per- 
formed at Medicina. Because blazars tend to be more variable at 
higher radio frequencies and a new 43 GHz receiver at Noto be- 
came available, from early 2006 data are acquired at 43 GHz in 
Noto (Bach et al. in prep.), while the radio monitoring at 22 GHz 
is done with the Medicina antenna only. 

The nominal system equivalent flux densities (SEFD) of the 
Medicina antenna are 296 Jy, 284 Jy, and 13 18 Jy at 5, 8, and 
22 GHz respectively, and 260 Jy, 770 Jy, and 800 Jy at the Noto 
antenna. On both telescopes the data is taken from the VLBI 
backend. Because of radio interference (RFI) the bandwidth at 
5 GHz is limited to 32 MHz at Medicina, but at the other frequen- 
cies the full IF (intermediate frequency channel) of 400 MHz 
bandwidth is used. 

The observations in Medicina are performed using ON/OFF 
measurements through the VLBI field system. The ON/OFFs are 
obtained by pointing the telescope subsequently ON the source 
and five beam widths to the east and west OFF the source. A typ- 
ical ON/OFF sequence (8 to 20 ON/OFF subscans) takes about 
8 to 18 minutes, depending on the frequency and source flux, 
and corresponds to 40 to 100 seconds ON-source integration 
time. Regular sky dip measurements (every 3 to 4 hours) are 
used to correct the data for opacity effects. The 22 and 43 GHz 
flux measurements in Noto are obtained with a similar method 
but the ON/OFF subscans are realized by "chopping" the sec- 
ondary reflector between the ON and OFF position, which al- 
lows a faster switching. Here, a typical ON/OFF sequence of 20 
subscans with 1 second integration time takes about 3 minutes. 
At the lower frequencies (5 and 8 GHz) the antenna performs 
raster scans in right ascension and flux densities are obtained by 
Gaussian fitting of the right ascension profiles. Opacity measure- 
ments are taken every half an hour. 

Since the VLBI field system is used to perform the observa- 
tions in Medicina and Noto, SNAP-command scripts or sched- 
ules are used to control the telescope. Most of the observations 
at Medicina are carried out remotely from Torino. Usually a sin- 
gle schedule is used for each 2-3 day observing run, switching 
between frequencies every few hours, but in case of technical 
problems or dramatic weather changes it is possible to change 
the schedule at any time. All the information provided by the 
observing system (e.g.: weather information, system tempera- 
ture, as well as the counts from the ON/OFF measurements) 
are written in ASCI format into a log-file. The field system also 
performs an analysis of the ON/OFF scans and provides uncal- 
ibrated antenna temperatures. However, to be more flexible in 
the data editing and to discard also single subscans, we wrote a 
Python program which extracts the information from the log-file 
and analyses again the individual subscans. The resulting table 
contains Stoke I antenna temperatures averaged over the two IFs 
containing the left and right circular polarized signals. The soft- 
ware which was used during a previous blazar radio monitor- 
ing pr oject performed in Medicina and Noto between 1996 and 
2000 dVenturi et al.ll200lb was not usable anymore, because the 
ON/OFF procedure of the field system has changed in the mean- 
time. Noto measurements from both observing methods are writ- 
ten in separate log-files and are reduced using IDL programs. 

The flux density measurements are calibrated using a 
set of standard primary calibrators, namely 3C 123, 3C 147, 
3C286, 3C295, 3C48, and NGC7027, and in addition at 
22 GHz DR21 and W30H. Their flux de nsities were calcu- 
lated by the po l ynom ial expressions given by lBaars et al ] (11977b 
and lOtt et all (Il994l) . Presumably the most stable calibrator 
is 3C286, for which 7.32 Jy (5 GHz), 5.22 Jy (8 GHz), and 
2.50 Jy (22 GHz) were adopted for the calibration, but at 22 GHz 



mainly NGC7027 (5.58Jy) and DR21 (19.04 Jy) were used, 
because they are more reliable (less variable) than quasars at 
this frequency. The calibrators are also used to obtain elevation- 
dependent gain curves to correct the data for this effect. The 
accuracy of the ON/OFF measurements in Medicina at 5 and 
8 GHz is usually about 5% to 10%, depending on the source 
flux density and weather conditions. At 22 GHz the accuracy 
is more weather-dependent, and ranges between 10% to 15% 
during good weather, but can be much worse during unstable 
weather conditions. Raster scans, like they are used in Noto, 
yield a slightly better accuracy, the typical error at 5 and 8 GHz 
being 2-5%. At 22 GHz the typical error at the Noto antenna is 
about 10%. 

3.2. Optical and near-IR data 

Contemporaneously to the radio observations we organise opti- 
cal {BVRI bands) and near-IR {JHK bands) observations at sev- 
eral observatories, which are listed in Table [2] The optical and 
near-IR magnitudes are obtained by differential photometry to 
known reference starfl 

The images were reduced using standard procedures for pho- 
tometry. To avoid offsets due to different calibration choices, the 
measurements from the different observatories were collected as 
instrumental magnitudes of the source and references stars and 
were calibrated all together. The uncertainty was calculated from 
the scattering of the references stars, which are supposed to be 
not variable. 



Table 2. The participating observatories. Given are the name, 
size of the telescope, observing bands, and number of observa- 
tions. Note that the data of Crimea and St. Petersburg are counted 
together, since they were observed by the same working group 
and were reduced all together using the same procedures. The 
uppercase numbers at the observatories denote the correspond- 
ing institutes (see page[T]l. 



Observatory 


Size [m] 




Bands 


N 


Medicina 12 


32.0 


5, 


8, 22 GHz 


2953 


Noto 5 


32.0 


5. 


8, 22 GHz 


1405 


Maidanak (T-60) 10 


0.6 




R 


3327 


Abastumani" 


0.7 




BR 


342 


Crimea 4 


0.7 




BVRI 


791 


St. Petersburg 


0.4 




BVRI 


Perugia 2 


0.4 




VRI 


80 


Campo Imperatore (AZT-24) 9 


1.1 




JHK 


2898 


Torino (OTAP) 1 


1.0 




BVRI 


1726 


Sabadell 7 


0.5 




R 


204 



Flux densities were derived from magn itudes using the ab- 
solute calibrations of iBessell et al] (1 19981). after d eredd ening 
according to the extinction laws of ICardelli et all ([1989). We 
adopted the B-band values of Galactic extinctio n provided by 
Nasa/I PAC extragalactic database^ (NED), from Schlege fet al] 
(Il998h . Non standard procedures that were used to reduce the 
data of 0235+164 and 2200+420 (BLLac) are described in 
Sect.|5] 



1 The adopted photometry can be found at 
http : //www. to . astro . it/blazars/webt/bamp/list .html 

2 http://nedwww.ipac.caltech.edu 



4. Results 

In this section we present the radio as well as the near-IR and op- 
tical light curves of the best sampled sources. The correspond- 
ing SEDs constructed with simultaneous data are also shown. 
Peculiarities in the light curves of single objects are described in 
the next section (Sect. [5]l 

4. 1 . Radio light curves 

Radio light curves at 5, 8, and 22 GHz of those blazars that are 
monitored from the beginning of our project are shown in Fig.Q] 
As visible in the legend, different symbols refer to different fre- 
quencies, whereas open symbols denote data taken in Medicina 
and filled symbols denote data from the antenna in Noto. The 
very good agreement between the radio measurements from the 
antennas in Noto and Medicina, taken with different methods 
and reduced by different authors, confirms a overall good qual- 
ity of the data. During some epochs the data was affected by bad 
weather conditions, which results in larger scattering and error 
bars, specially at 22 GHz. 

Most of the sources display significant variability on 
monthly time scales. Some sources show an achromatic be- 
haviour, where the flux densities at different frequencies vary si- 
multaneously, e.g. 1510-089. In other cases variations at 22 GHz 
are likely leading those at the lower frequencies, e.g. 1226+023 
(3C273) and 2251 + 158 (3C 454.3). In the case of 2200+420 
(BL Lacertae), we recognise both characteristics. A flare that is 
more pronounced at 22 GHz appeared in early 2005 and a second 
one equally bright at all frequencies followed about 250 days 
later. 



4.2. Optical light curves 

Optical light curves are not shown for all sources, but only for 
those which have a good time coverage (Fig. 01. The lack of data 
is mostly due to solar conjunction or bad weather periods. The 
measurements from different observatories agree quite well and 
the variations in the optical BVRI bands and near-IR JHK bands 
usually show similar amplitudes. To allow the reader to better 
distinguish between different bands only B, R, J, and K light 
curves are shown. Measurements in V, I, and H are shown in the 
SEDs only (next Sect.). In general, the sources are all more vari- 
able in the optical bands than at radio wavelengths. In Fig. [2] we 
complement the light curves of 0235+164 and 2251 + 158 with 
data taken by the Whole Earth Blazar Teles cop (WEBT) dur- 
ing specific campa igns (grey symbo ls, see iRaiteri et al.l [2006b 
for 0235+164 and lVillata et aDl2006l for 2251+258). 

Part of the data presented here were also provided to other 
multi-wavelength campaigns and the specific time intervals are 
marked by yellow (grey in the printed version) vertical bars in 
the figure. For 0851+202 the periods of two short-term WEBT 
campaigns organized in conjunction with two pointings of the 
X-ray satellite XMM-Newton (Ciprini et al. in prep.) are high- 
lighted. Similarly, the yellow strip crossing the 1253-055 light 
curves indicates the time of the January-March 2006 WEBT 
campaign on this object (Bottcher et al. in prep.). In the case 
of 2251 + 158, the yellow area corresponds to the new WEBT 
campaigns following the 2004-2005 one (Villata et al. in prep; 
Raiteri et al. in prep.). 



3 http : //www . to . astro . it/blazars/webt/ 



4.3. Spectral energy distributions 

In Fig.[3]the SEDs of several blazars are shown. Grey dots rep- 
resent non-simultaneous literature data taken from NED, while 
different symbols (and colours) have been used to distinguish 
the contemporaneous data acquired for this project. The sources 
shown here are those for which several epochs of simultaneous 
data exist and whose SEDs show significant variability. In this 
energy range, we expect to see the synchrotron component peak- 
ing in the sub-millimetre to infrared region, in the case of these 
low-energy peaked blazars. 

The lines connecting the data points are third-order polyno- 
mial fits, and do not represent any emission model. Although 
the SEDs of blazars seem to be represented quite well by log- 
parabolic fits, as demonstrated by various aut hors dLandau et al.l 
1986; iPerri et all [20031; iMassaro et alll2004allbh . we used third- 
order fits because they show more clearly the optical spectral 
changes which are present in some sources. In those cases where 
a "blue bump" appears, we used cubic spline interpolations. 
Thus the lines are drawn just to guide the eye and to better il- 
lustrate the behaviour of the source. A proper modelling of the 
SED is behind the scope of this article, but will be presented in 
a future publication. Therefore, the peaks of the SEDs in Fig. [3] 
might not exactly match the position where most of the energy 
is released, but should provide an estimate of the true location. 

As already seen in the radio light curves, also the broad-band 
variability can be either chromatic or achromatic. The long-term 
variability of the BLLac object 0716+714 seems to affect the 
emission from the radio to the optical wavelengths at the same 
time, whereas the behaviour of 2251 + 158 is completely differ- 
ent. Indeed, in 2005 the source experienc ed a historical flare , 
peaking around May in the optical b and dFuhrmann etaD[2006; 
iPian et al.ir2006b IVillata et aLll2006l) . and possibly a bit later in 
the near-IR one, while it appeared at 22 GHz with a few months 
of delay, in August-September 2005 (see Fig.HJ. 

Although most blazars of our sample reach their emission 
peak in the millimetre-to-infrared regime and their SED declines 
in the optical band, some sources show different behaviours. 
Exceptions are e.g. 0716+714, whose SED appears rather flat 
in the optical band, and the quasars 0836+710, 1226+026, and 
1510-089, whose optical part of the SED is increasing again 
(see the individual paragraphs in Sect. [5] for more details). 
Commonly, this increase is referred to as a "blue bum p", which 
could be due to emission from an accretion disk (e.g.lBregman 
Il990h . 

For some sources, e.g. 0235+164, 1156+195, and 
2251 + 158, we can also see spectral changes during differ- 
ent states in the optical, whose dependence on the brightness 
level, if any, is not clear. In the other cases the colours do not 
change much during different brightness levels. 

5. Notes on individual sources 

0219+428 (3C 66A) This BL Lac obj ect has been observe d in- 
tensively at optical wavelengths (e.g. iBottcher et al.ll2005l) and 
lLainelaet al.l (il999) found a possible periodicity of 65 days in 
its optical bright stat e. Very fast superlum inal motion of up to 
29 c was reported by Jorsta det al.l d200ll) . Despite some small 
flux changes, the radio light curves reveal not much variability 
(Fig. [J}. A major problem for radio measurements of 3C 66A is 
its close companion 3C 66B, which at some position angles ap- 
pears in the OFF measurements or the cross scans and makes the 
calibration more difficult. During the first ~ 50 days of the op- 
tical light curve a brightness decrease is visible, which reached 



its minimum at about JD=2453650 (October 2005), with some 
shorter variability on top of the long-term trend (Fig.©. 



0235+164 This BLLac object usually displays strong variabil- 
ity at optical and radio wavelength. Although very compact, 
sometime s jet features with high superluminal motion can be 
observed. lJorstad et al.l ((2001) reported an apparent speed as fast 
as 40 c, and anot her feature at about 26 c was observed by 
Pine r et al. (2006). An outburst in early 2004 was expected as 
part of a ~ 5.7 year quasi-periodicity resulting fro m the analy- 
sis of the radio (and optical) historical light curves (iRaiteri et al.l 
120011) . A long-term multi-wavelength campaign wa s organized 
by the WEB T to monitor the predicted outburst dRaiteri et al.l 
2005, 12006allbh . The grey symbols in Fig. [2] represent the data 
taken in the last period of the WEBT campaign. Only two arc- 
seconds south of the source, there is another AGN, which signif- 
icantly affects the optical flux measurements of the 0235+164. 
Therefore, the magnitudes shown in Fig.|2]have been corrected 
for this contribution. Moreover, the near-IR and optical fluxes 
shown in Fig. [3] were obtained by taking into account not only 
the Galactic extinction, but a lso the extinction d ue to the fore- 
ground galaxy (see details in IRaiteri et al J 12005). We observed 
the source in a very low state. The light curves show a slow 
increase at radio frequencies, reaching a maximum at 22 GHz 
around JD=2453700 (November 2005) and several weeks later 
at 8 and 5 GHz. Some significant variability is also visible at the 
optical wavelengths. 



0716+714 This BLLac object is one of the best-studied intra- 
day variable sources. Simultaneous radio, optical, UV, and X-ray 
monitoring yielded a short duty cycle of variability at all fre- 
quencies, as well as a correlat ion between the r apid variations 
at different frequency regimes (Wagner et al. 1996). Those mea- 
surements suggested very small sources sizes and high bright- 
ness temperatures exceeding 10 17 K. The data collected dur- 
ing a recent large multi-wavelength campaign allowed to ob- 
tain more moderate valu es of 10 14 to 10 16 K (lOstorero et al.1 
2006; lAgudo et al.l [2006; Fuhrmann et al. in prep.), well in 
agreement with the Dop pler factors derived from kinemati- 
cal studies of the jet (e .g, 1Jorstad et aT1l200lb iKellermann et al.l 
2004; lBachetal.ll2005l) and sp ace VLBI measurements of the 
all 



source size (Bach et all [2006 aj). A detailed study of the multi- 
wavele ngth long-term variability was performed by Raiter i~et~ai1 
(120031) . During the present monitoring, 0716+714 showed three 
flares in the radio, lasting about 100 days, and exhibiting a 
smaller amplitude than those occurred in the past. But in our 
observing period, the radio flux also went down to its histori- 
cal minimum of about 0.5 Jy (Fig. [TJ. Similar variations appear 
also in the optical band (Pig. [2j, and an inspection of the SEDs 
confirms that the optical and radio flux densities vary contem- 
poraneously (Fig. |3). The fact that the peak of the SEDs falls 
at higher energies compared to many of the other blazars in our 
sample confirms the classification of this source as an "inter- 
mediate" BL Lac (as 3C 66A), half-way between the low-energy 
peaked (LBL) and the high-energy peaked (HBL) BLLacs. 



0836+710 This ultra-luminous quasar at z — 2.17 is show- 
ing a well correlated broad-band variability that also seems to 
correlate with the st ructural changes seen in the VLBI jet (e.g. 
lOtterbein et al.llT998h . The source showed only weak variability 
(Figs. Q] & |2]i. The hard near-IR-optical spectrum implies that 
the blue bump component, likely due to thermal emission from 



the accretion disc, is still dominant in the near-IR band, and con- 
sequently it is not possible to distinguish the high-energy part of 
the synchrotron component. This is actually expected, because 
of the high redshift of the source. 



0851+202 (OJ287) This well-studied BLLac object is known 
for its ~ 12-year period icity in the optical light curve 
(|Sillanpaa et all 1 1 988l 1 1 996l) . The outbursts appear as double- 
peaked and their periodic occurrence is mos tly interpreted 
in terms of a binary black hole system (Sillanpaa et al 
1988: lLehto & Valtonenlll996HVillata et alJl99SlValtaoia et al 
2000). Depending on the model, different dates around 2006 



are predicted for the next outburst, but the rece nt strong activ- 
ity mak es an identification not easy. According to lValtonen et al.l 
(2006a b), the outburst occurred in late 2005 may have already 
been the predicted one, but a final statement is not possible be- 
fore 2008, when all predicted dates will have passed. During our 
radio monitoring we observed an increasing trend, with some 
flares detectable at 22 GHz (Fig. [TJ. The optical light curves 
show the declining phase of the large outburst around November 
2005, and another flare around April-May 2006 (Fig. 0. The 
SEDs shows that the amplitude of the optical variations of 
OJ 287 were slightly larger than in the radio (Fig. [3j. 



0954+658 Th i s BL L ac object is another fam o us IDV source 
(IWagner et al.l 119931; IWagner & Witzell 119951; iGabuzda et all 
2000). The jet structure shows many bends from milliarcsecond 
to arcs econd scales and speeds of up to 9 c were measured with 
VLBI ( IGabuzda & Cawthornell996tlKollgaard et al.lll992l) . The 
long-term radio light curve (Fig. [TJ displays some variability, 
which is more pronounced at 22 GHz, while greater activity 
characterizes the optical behaviour (Fig. |2j. This becomes even 
more evident in the SEDs, where in addition one can notice the 
flattening of the optical spectrum during the flare of January 
2005 (Fig. [3]l. 



1156+295 This quasar is extremely variable over al l the elec- 
tromagnetic spectrum, from radi o waves to y-rays dWills et al.l 
119831 ll992trHartman et al.lll999l) . 1156+295 has a straight arc- 
second scale jet and a more curved jet on milliarcsecond scales, 
with j et compone nts moving at superluminal speeds of up to 16 c 
(e.g. Hong et al .1120041) . Our observations reveal a major optical 
outburst in March 2006 (Fig. [2}. An increasing trend is also vis- 
ible in the radio light (Fig. [TJ and it will be interesting to see if 
we can find a connection between the two events. 



1226+023 (3C273) This object shows all the characteristics 
that are typical of high-luminosity quasars: a flat radio spec- 
trum of the core, str ong and rapid variability in all the observed 
energ y ranges (e.g. Ivon Monti gn y et al.l 119971: iMcHardy et al.l 
1999), variabl e polarization, and a radio jet with superluminal 
motion (e.g. lPearson et alJll98ll) . Additionally, it shows an op- 
tical and X-ray jet, and a very prominent UV excess, the so - 
called big blue bump (e . g. iRoser et al.ll 1 996t I Jester et al JI2005I) . 
In the SEDs the contribution of this component begins to domi- 
nate over the synchrotron one in the optical band (Fig. [3j. Since 
it is likely due to thermal radiation from the accretion disc, it is 
not expected to be strongly variable on short time scales, and in- 
deed the optical light curves do not show important flux changes 
over the whole time range (Fig.[2]i. On the contrary, a noticeable 
flux increase has recently been observed at 22 GHz, which might 



propagate also to the lower radio frequencies in the next future 

(Kg.©. 



the SEDs (Fig. fJJ, specially when considering that the last four 
SEDs are separated by only a few weeks. 



1253-055 (3C279 ) This quasar, c onsidered the first superlu- 
minal source dWhitney et al.l Il97ll) . is often nearly as bright 
as 3C2 73 at high frequenci es. Since it is also a bright y-ray 
source dHartman et al.l ll999). 3C279 has been the target of a 
number of campaigns of cont emporaneous monitoring from y- 
rayto radio frequencies (e.g. iGrandi et al.l 11996k IWehrle et al.1 
1998). Another campaign has recently been performed, which 
was motivated by an optical bright state observed in December 
2005-January 2006, which in turn triggered ToO observations 
by Chandra and INTEGRAL (Bottcher et al. in prep.). The ra- 
dio light curve at 22 GHz shows a noticeable decrease in the first 
months of the project, and a more moderate increase afterwards; 
on the contrary, the radio flux at 8 and 5 GHz did not change 
very much (Fig. [TJ. Due to the low declination of the source, 
the optical light curve has a large gap around the solar conjunc- 
tion (Fig. 0. The SEDs show a fast changing optical spectrum 
(Fig.0. 



1510-089 This is another extreme quasar, highly polarized, and 
bright at X-ray and y-ray energies. It exhibits a highly bent 
radio jet structure with fast superluminal speeds of up to 20 c 
dSingh et al.lll997t iHartman et al.lll999b Soman et al.ll200Tl) . In 
the radio band, we saw the source fading simultaneously at all 
frequencies (Fig. [TJ. Noticeable variations were seen also in the 
near-IR (Fig. |2). Note also that our SEDs (Fig. [3) show a dip 
between the ne ar-IR and optical bands, in agreement with the 
observations by lNeugebauer et all dl979l ) in 1977 (data reported 
by NED, grey dots). In that case, however, the optical spectrum 
was quite harder and it seems that a much stronger emission 
component was present towards smaller wavelength. Such an in- 
crease in the spectrum of a quasar is comm only referred to ther- 
mal emission from an accretion disk (e.g. lBregmanl 19901) . but in 
that case, although the observations are separated by ~ 30 years, 
the component should not be so strongly variable. A similar be- 
haviour was observed in AO 235+164, where a variable UV- 
to-soft-X-ray bump was found dRaiteri et al.ll2005l l2006allbh . In 
this case the authors speculate that the variability could be due 
to a strong change of e.g. the accretion rate. Another possibil- 
ity would be the presence of another emission component, e.g. 
synchrotron radiation from an inner jet region with respect to 
that producing the radio-near-IR radiation. For 1510-089 the 
latter explanation would be favoured by the notic eable optical 
variab ility also observed on short time scales (e.g. lRaiteri et al.l 
1998). Further multi-wavelength monitoring will help to better 
characterize the behaviour of this component. 



1641+399 (3C 345) The quasar 3C 345 has been studied exten- 
sively at all accessible wavelength. The observed strong vari- 
ability in the optical and radio bands, the existence of possible 
periodicities, and the detection of bent trajectories of the radio 
jet components, made this source particularly inter esting to test 
jet pr e cession and binary black ho l e models (e.g. iBiretta et al.l 
1986; Caproni & Abraham 200 41 lLobanov & Roland! 120051: 
Lobanov & Zensusl 1999; Schramm et al.l 119931: IStevens et al.l 



2200+420 (BL Lacertae) This is the prototype of the BL Lac 
objects, which are well-known for their pronounced variability 
at all wavelengths, from the radio to the y-ray band. BLLac 
has been studied intensively since its discov ery and was the tar- 



( Bloom et al. 



Villata et al.1 



1997 



2002 



1996; Webb et al.l ll994l Our radio light curves are slowly de- 
clining (Fig. |TJ, but in the optical band a strong increase can be 
seen from the beginning of the last observing season (Fig. 
This becomes even more evident by looking at the evolution of 



get of several m ulti-w avelengths campaigns 
Samb runa et alJI 19991 iMadeiski et alj|l999t 
iBottcher et al.l 120031; IVillata et al.ll2004allbl) . The radio emission 
appears well correlated, where variations at higher frequencies 
lead the lower- frequency ones by several days to a few months 
fe.g. IVillata" et al.ll2004al) . The detection of a fair correlation be- 
tween the optical variations and the radio variability from the 
VLBI- core suggests a co mmon origin in the inner portion of 
the jet dBach et al.ll2006bl) . A number of superluminal jet com- 
ponents have been observ ed displaying bent trajectories and 
speeds from 3 c to 9 c (e.g. iMutel et aljfl990t iDenn et al 12 000: 
IStirling et al]l2003l iKeUermann et al.ll2004t iJorstad et alJl2005h . 
During our monitoring BL Lac showed its typical variability be- 
haviour (Figs. [JJ & [2j. There are flares which are more pro- 
nounced at higher frequencies (including optical counter parts) 
and appear delayed at lower frequencies, and flares that are oc- 
curring contemporaneously and with comparable amplitudes, 
but appear only at radio wavelength. To construct the SED the 
optical flux densities of BL Lac were corrected for the contribu- 
tion o f the host galaxy using the /?-band value by Scarp aeTail 
d2000l) (R = 15.55 + 0.02) and deriving the other bands from 
the colour indices for elliptical galaxies with My < -21 from 
iMannucci et al.l d200ll) . 



2230+114 (CTA 102) This quasar shows the typical blazar be- 
haviour and appears in the source lists of many radio and opti- 
cal monitoring projects. Multi-wavelength monitoring and imag- 
ing revealed correlated variability and radio jet speeds of up to 
21c w ere derived from VLBI observations (e.g. Rantakvr o et al.l 
120031) . Our near-IR-optical light curve is dominated by a huge 
flare lasting more than a hundred days around October 2005 
(Fig.|2]i, and we can follow now an ongoing radio flare (Fig. [TJ. 
The SEDs show impressively the change of the spectrum with 
the propagation of the flare towards lower frequencies (Fig. fJJl, 
similarly to the case of 2251 + 158 (see next paragraph). 



2251+158 (3C 454.3) This is another well-known quasar and 
also one of the brightest ones. The source is showing re- 
markable variability at all wavelengths, high degrees of po- 
larization, superl uminal motion of t he radio jet components , 
and a X-ray jet dCotton et al.lll984t iPauliny-Toth et al.1 119871: 
IHartman et al.lll993l) . With the start of our monitoring 3C 454.3 
underwent an exceptional optical outburst lasting for more than 
one year, subsequently fol l owed by a radio fl are (Figs. [JJ & El 
and iFuhrmannetai] 120061; IVillata et al.l 120061) . The maximum 
brightness detected was R = 12.0, which, at z — 0.859, rep- 
resents the mo st luminous quasa r state observed until today 

(M B 31.4. IVillata etai] 120061) . Therefore, the SED shows 

a remarkable variability and the propagation of the flare from 
higher to lower energies is clearly visible (Fig.fJJl. At the begin- 
ning of the rising phase the near-IR-optical spectrum appears 
rather flat and changes to a very steep one after the brightness 
peak was reached. The radio flare did not reach frequencies 
much below 15 GHz and, although some structural changes in 
the radio jet were detected, so far no jet component related to 



the outburst is visible (Marscher et al. 2007, in prep). A detailed 
analysis of the radio-optical correlations will be presented in an 
upcoming paper (Villata et al. in prep.). 

6. Conclusions 

During the first 18 months of our multi-wavelength monitoring 
we obtained about 4000 data points in the radio bands, at 5, 8, 
and 22 GHz, and 5500 near-IR and optical measurements. 

Some of our sources were targets of observations by X- 
ray satellites and by ground-based TeV telescopes, and we 
contributed data to the multi-w avelength WEBT campaigns on 
0235+164 dRaiteri et al.ll2006bl) . OJ287 (Ciprini et al. in prep.), 
3C279 (Bottc her et al. in prep.), BL Lac (Villata et al. in prep.), 
and 3C 454.3 (IVillata et al]l2006t Villata et al. in prep; Raiteri et 
al. in prep.). 

Most of the sources display significant variability on 
monthly time scales. Some sources show an achromatic vari- 
ability behaviour, where the frequencies vary almost simultane- 
ously, as e.g. 1510-089. In other cases the flux changes at the 
higher frequencies are likely leading lower-frequency ones, as 
for 1226+023 (3C273) and 2251 + 158 (3C 454.3), where flares 
are visible at 22 GHz, but are not yet visible at 8 and 5 GHz. In 
the case of 2200+420 (BLLac) we find both characteristics. A 
flare that is more pronounced at 22 GHz appeared in early 2005, 
and a second one equally bright at all frequencies appeared about 
250 days later. 

Many possible mechanisms which are able to explain the 
various variability behaviours are discussed in the literature (see 
Sect. [T|for references). They manly divide into geometrical mod- 
els where changes of the viewing angle lead to variations of the 
Doppler beaming factor and models which propose changes of 
the physical parameters of the jet plasma itself. Because of the 
limited time interval covered by the present project, it is not pos- 
sible yet to favour or rule out some of these models. Neither does 
it allow to perform a meaningful statistical analysis on the time 
scales of variability or on cross correlations among flux changes 
at different frequencies. However, the planned continuation of 
this multi-wavelength monitoring during the next years is ex- 
pected to provide a wide database for our targets that will allow 
us to extract more information on the mechanisms at the origin 
of blazar emission variability. 

Moreover, with the launch of AGILE and GLAST the GeV 
energy domain will be accessible and our monitoring will sup- 
port the high-energy observations with low-energy data, to learn 
more about the inter-connection between the synchrotron and 
inverse-Compton components of blazar SEDs. 
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Fig. 1. Radio light curves at 5, 8, and 22 GHz of our blazar sample between December 2004 and June 2006 taken with the antennas 
in Medicina (open symbols) and Noto (filled symbols). Although some of the 22 GHz measurements suffer from bad weather 
conditions, we are able to trace the variability behaviour very nicely. In most cases the flares appear first at the higher frequency, 
suggesting that the variability is caused by a shock propagating along the jet (continued on next page). 



2005 2005.5 2006 2005 2005.5 2006 



4 — 



3 - 



3 



I — 



1/3 



20 
18 
16 
14 
12 

5.5 
5 

4.5 
4 



n 1 ' | 

0851+202 



J I I I I I I I L I L 



1156+295 



% 

in 



CD [f 



4>" 



ti — — in 



— - ♦' 



J I I I I I I I I I L 



1253-055 



•**J 



J I I I I I I I I I L 



1611+343 
oil 



f 



I 



<JH> 



II'. 



* 



J I I I I I I I I I L 



"1 1 1 — 

0954+658 



— 2 



♦ \ 



-.El-y-f 



•Dcs' Vnii 



¥ 1 



— 0.5 



J I I I I I I I I I L 



1226+0^3 



■<>■#■■■«))" 



•f f 



J I I I I I I I I I L 



1510-089 



J I I I I I I I I I L 



1633+382 



- k \M ^ * *. 

_i I I I I I I I I I L 



o 5 GHz 
□ 8 GHz 
1-5 o 22 GHz 



36 
32 
28 
24 
20 
16 



on 



— 2 



— 4 



8400 8500 8600 8700 8800 8400 8500 8600 8700 8800 

JD-2445000 



Fig. 1. continued. 
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Fig. 2. Near-IR ( JK) and optical (BR) light curves of those source where we could get a good time coverage. The different colours 
(only in the electronic version) and symbols denote the different observing bands. Grey symbols denote already published data and 
yellow strips indicate the periods when WEBT campaigns were active, see Sect l4.2l for more details. Note that the time range shown 
is not the same for all sources. 
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Fig. 2. continued. 
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Fig. 3. Examples of blazar SEDs, which we obtain from our data for those sources which are best sampled. Grey dots represent 
archival data taken from NED. See Sect.[5]for more details. 
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Fig. 3. continued. 



